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NOVEL PHARMACOPHORE FOR THE DISCOVERY AND TESTING OF NA, 
K-ATPASE INHIBITOR COMPOSITIONS AND METHODS FOR THEIR USE 
IN TREATING CARDIOVASCULAR DISEASES AND CONDITIONS 

CROSS REFERENCE TO RELATED APPLICATIONS 
[0001] This application is a continuation-in-part of United States Provisional Application 
Serial Nnmber 60/425,037, filed November 7, 2002. 

FIELD OF THE INVENTION 
[0002] The present invention relates to the discovery of a three dimensional 
pharmacophore based on a model of the Na, K-ATPase that will enable the discovery and 
design of novel positive inotropic agents using molecular modeling, pharmacophores, rational 
dmg design techniques, and methods of using the novel compoimds in treating heart diseases. 

BACKGROUND OF THE INVENTION 
[0003] Each of the references cited herein is incorporated by reference in its entirety. 
[0004] Heart diseases, such as cardiac failure, paroxysmal atrial tachycardia, and 
caridomyopathy, are a major cause of morbidity and mortality aroxmd the world. Cardiac 
failure occurs when the heart becomes unable to pump blood effectively at a rate that meets 
tiie needs of the metabolizing tissues. The underlying cause of cardiac failure lies in the 
reduced contractility of heart muscles. Reduced contraction of the heart leads to reduced 
heart output of blood volume while the amount of blood being retumed to the heart remains 
the same. The residt is an increase in heart blood volume, which in turn leads to heart 
congestion. Congestive heart failure causes lowered blood pressure, poor renal blood flow 
with concomitant edema in the lower extremities and pulmonary edema, as well as renal 
failure. Positive inotropic activity, meaning increasing the force of the contraction of the 
heart, is thus a primary goal in the treatment of heart failure patients. 
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[0005] While recent advances in diagnostic and surgical procedures have revolutionized 
cardiovascular medicine, similar progress has not been attained with respect to 
pharmacological agents. Increasing the force of contraction of the failing heart is still largely 
achieved by cardioglycosides (also called cardioactive glycosides) derived from the plant 
Digitalis purpurea, (foxglove) such as digoxin and digitoxin. Unfortunately, digoxin, 
digitoxin, and other cardioglycoside drugs, the most widely used clinical inotropic agents, 
have narrow therapeutic indices and a high potential for considerable toxic effects mcluding 
ventricular arrhythmias. Therefore, the difference between the dosage leading to a therapeutic 
increase in heart muscle contraction and the dosage leading to fatal toxicity is minute. This 
potentially deadly property severely curtails the clinical usefidness of the currently available 
dmgs. While much effort has been directed toward the discovery of novel therapeutic agents, 
the lack of three-dimensional stmctural coordinates for the receptor of inotropic drags. 
Sodium Potassium (Na, K) ATPase, has markedly limited success. Increasing the therapeutic 
index of positive inotropic Na, K-ATPase (also referred to as Na+, K+-ATPase hra-ein) 
inhibitors would thus represent a major advance in the clinical management of heart disease. 
[0006] Sodium Potassium ATPase is a heterodimeric transmembrane protein that actively 
exchanges sodium and potassium ions across cell membranes using energy from the 
hydrolysis of ATP. Inhibition of the Na, K-ATPase in the myocardial cells of the heart 
muscle results in the sequential increase of cytosolic sodium and calcium concentrations, 
thereby increasing the force of contraction. Regulation of the levels of sodium and potassium 
contribute to many essential physiological processes including maintenance of the membrane 
potential for muscle contraction. 

[0007] Cardioglycoside inhibitors bind to the extracellxilar domain of the catalytic 1 a 
subunit of the Na, K-ATPase and inhibition of the ATPase results in both the beneficial and 
detrimental effects associated with cardioglycoside administration. While researchers have 
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obtained low resolution crystal structures (--^A) of the Na, K-ATPase and gained structural 
insights from solid-state NMR studies of ATPase-inhibitor binding, the dearth of detailed 
structural information has hindered flie discovery of novel inotropic inhibitors. One of the 
main obstacles to Na, K-ATPase crystallization is that the protein is structurally complex, 
possessing ten transmembrane spanning domains (Hl-HlO), which makes the development of 
the high quality crystals of integral membrane proteins required for X-ray crystallography 
exceedingly difficult. 

[0008] The exact mechanism of inhibition has remained elusive. Although researchers 
have examined altemative inhibitory agents, of the hundreds of cardioglycosides tested, none 
have been shown to have larger therapeutic indices than existing therapeutic agents do. 
Moreover the lack of structural knowledge of the Na, K-ATPase has thwarted the structural 
investigation and structure based discovery of novel altemative treatments. A novel ^proach 
to the design, discovery, and use of inotropic dmgs with reduced toxicity would thus be very 
useful in the treatment of heart failure. 



SUMMARY OF THE INVENTION 
[0009] The present invention creates a novel pharmacophore from which to generate 
scaffolds, and then molecules, which have inotropic properties. Preferably, the 
pharaiacophore is defined by the parameters of Table 4 and Table 5, which reference the 
spheres in Figure 4. The invention also contemplates the method of creating the novel 
pharmacophore. 

[0010] Another aspect of the present invention sets forth a basic scaffold created using 
the novel pharmacophore. The basic scaffold molecule is a novel Na, K-ATPase Ugand. The 
scaffold can be used to develop positive inotropic drugs that inhibit Na, K-ATPase activity. 
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[0011] A third aspect of the present invention discloses further developed potential 
scajSblds of the novel Na, K-ATPase ligands. These scaffolds represent the parent molecules, 
which satisfy the basic three point pharmacophores. The Na, K-ATPase inhibitor compounds 
are of the formula: 




wherein Rl, R2, R3 and R4 can be any organic functional group containing a hydrogen bond 
donor or a hydrogen bond acceptor and X is any element or group that allows the compound 
to retain inotropic activity. Preferably, X is N, O, S, or C. 

[0012] A further aspect of the present invention takes the scaffold molecules and 
develops them into fimctional dmgs for Na, K-ATPase inhibition, which can be used to 
prevent or treat heart conditions, particularly congestive heart failure and supraventricular 
arrhythmia. Eight compounds have shown promise as inotropic agents, with three of the 
eight demonstrating high therapeutic potential based on their ATPase inhibitory 
characteristics. 

[0013] These and other aspects of the invention will be further elucidated in the following 
descriptions and examples. The descriptions in the present invention are provided only as 
examples and should not be understood to be limiting on the claims. Based on the 
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description, a person of ordinary skill in the art may make modifications and changes to the 
preferred embodiments, which does not depart from the scope of the present invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0014] Figure 1. This figure shows how the threaded alignment of the la-subunit of 
SERCAla (skeletal muscle sarcoplasmic reticulum/endoplasmic reticulum Ca2'*"-ATPase) 
[SEQ ID NO: 1] with the sheep Na, K-ATPase [SEQ ID NO: 2] was accompUshed. The 
alignment shown is the result of "threading" the Na^, K^-ATPase sequence to that of the 
template. Identical residues between the two sequences are illustrated with filled black 
boxes, while similar residues correspond to the filled gray boxes. Regions modeled are noted 
by Unes above lettered designations, with the letter shown at the beginning and end of the 
sequence encoding the region. The letters correspond to letters attached to the helices and 
loops depicted in Figures 2 A and 2B, The solid lines under the sequence of Figure 1 represent 
helical structures, while ttie broken lines are indication of loop regions. Local areas of high 
identity correlating with functional motifs are shown by open black boxes and described in 
the text. Positions of sequence changes between the sheep and rat Na***, K^-ATPases are 
denoted with an asterisk (*). 

[0015] Figure 2. This figure shows the resulting homology model of the extracellular 
loops aad the transmembrane domains regions of the sheep la-subunit of the Na, K-ATPase. 
Shown is a visual representation of the three-dimensional coordinates assigned to the 
homology model (2A), and a representation of the a-heUces from the extracellular domain 
(2B). Only alpha carbons are depicted. The numbering, showing the orientation of the 
helices, is consistent with the information obtained from low-resolution structures. 
Tryptophan and lysine residues are displayed (2A). The lettering corresponds to the amino 
acid sequence portions shown in Figure 1. 
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[0016] Figure 3. This figure shows the binding orientation of cardioglycosides to the la 
subunit of the Na, K-ATPase predicted by the methods of the invention. White numbers in 
the figure indicate transmembrane domains. The extracellular (3 A) and cross-section (3B) 
represent the consensus orientation for all nine docked cardioglycosides. Both Na^, K^- 
ATPase and ligands are depicted with Connelly surfaces. The figure depicts consensus 
orientation of 25 individual dockings of ouabain to the la-subunit of the sheep Na^, K^- 
ATPase (3C). Orientation of digoxin (3D, upper) and digoxigenin bisdigitoxose (3D, lower) 
illustrate the identical conformation regardless of the number of sugar moieties. 
[0017] Figure 4, This figure shows the pharmacophoric elements fitting the binding 
criteria for la-subunit Na, K-ATPase inhibition. Each pharmacophoric element as described 
herein is represented by a point encased in a sphere. The spheres (diameter = 2.5 A) represent 
the tolerance of location allowed during database searching. The pharmacophoric elements 
are combined to form three, four and five point pharmacophores representing increasing 
stringency of search criteria (Table 4). 

DETAILED DESCRIPTION OF THE INVENTION 

Definitions 

[0018] The present invention can best be understood in light of the folloAving 

definitions. 

[0019] *Thotropic" means affecting the force or energy of muscular contractions. 

[0020] "Glycosides" are any compound that contains a constituent sugar in which the 
hydroxyl group attached to the first caibon is substituted with an alcoholic, phenolic, or other 
group. 
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[0021] '^Rational drug design" means modeling the molecular structure of the target of a 
drug and then designing a dmg that will attack tiie target. 

[0022] 'Tharmacophore" is a model for developing one or more molecular scaffolds or 
structure used as the basis for drug development. 

[0023] "Scaffold" means a supporting framework, such as the basic backbone or structure 
of a molecule witii designated positions capable of containing a one or more molecules or 
functional groups. 

[0024] *T)ynamic programming" is a use of scoring matrices and gap penalties to produce 
optimal alignment 

[0025] 'TTnaltered" means not changed, substituted, deleted, or mutated. 
[0026] "Pharmaceutically effective dose" as used herein refers to the amount of an novel 
inotropic compound or novel compound composition comprising one or more of the novel 
compounds described herein that produces a desired therapeutic effect, such as treating the 
target disease. The precise amount of the pharmaceutically effective dose of a novel 
compound or novel compound composition that will yield the most effective results in terms 
of efficacy of treatment in a given subject will depend upon the activity, pharmacoldnetics, 
pharmacodynamics, and bioavailability of a particular inotropic compound, physiological 
condition of the subject, including age, gender, disease type and stage, general physical 
condition, responsiveness to a given dosage and type of medication, the nature of 
pharmaceutically acceptable carrier in a formulation, and a route of administration, among 
other potential factors. Those skilled in the clinical and pharmacological arts will be able to 
determine these factors through routine experimentation consisting of monitoring the subject 
and adjusting the dosage. Remington: The Science and Practice of Pharmacy (Gennaro ed. 
20* edition, WiUiams & WiUdns PA, USA) (2000). 
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[0027] As used herein, the 'target disease" may be any disease in which inhibiting Na, K- 
ATPase treats the disease. Without limitation, the target diseases are heart diseases, 
particularly congestive heart failure and supraventricular arrhythmia. While it is possible for 
a novel compound to be administered as a pure or substantially pure compound, it is 
preferable that the novel inotropic compound be administered as a composition m the form of 
pharmaceutical formulations or preparations suitable for a particular administration route. A 
novel compound composition comprises one or more novel inotropic compounds and a 
phaimaceutically acceptable carrier. 

[0028] The term "phaimaceutically acceptable carrier" as used herein means a 
pharmaceutically-acceptable material, composition or vehicle, such as a liquid or solid filler, 
diluent, excipient, solvent or encapsulating material, involved in carrying or transporting one 
or more novel inotropic compounds from one tissue, organ, or portion of the body, to another 
tissue, organ, or portion of the body. The novel inotropic compounds may be administered in 
Uposomes. Each component must be "pharmaceutically acceptable" by being compatible 
with the other ingredients of the formulation. It must also be suitable for use in contact with 
the tissue or organ of humans and animals, specifically the heart and circulatory system, 
without excessive toxicity, irritation, allergic response, immunogenecity, or other problems or 
compUcations, commensurate with a reasonable benefit/risk ratio. 

[0029] A "route of administration" for a novel compound or composition can be by any 
pathway known in the art, including without limitation, oral, enteral, nasal, topical, rectal, 
vaginal, aerosol, transmucosal, transdermal, ophthalmic, puhnonary, and/or parenteral 
administration. A parenteral administration refers to an administration route that typically 
relates to injection. Parenteral administration includes, but is not limited to, intravenous, 
intramuscular, intraarterial, intraathecal, intracapsular, infraorbital, intracardiac, intradermal, 
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intr^eritoneal, transtracheal, subcutaneous, subcuticular, intraarticular, subcapsular, 
subarachnoid, intraspinal, via intrastemal injection, and/or via infusion. 
[0030] 'Treatment' • of or * treating" a disease may mean preventing the disease by causing 
clinical symptoms not to develop, inhibiting the disease by stopping or reducing the 
symptoms, the development of the disease, and/or slowing the rate of development of the 
disease, reUeving the disease by causing a complete or partial regression of the disease, 
reducing the risk of developing the disease, or a combination thereof 

[0031] The term "contacted" when applied to a cell, tissue or organ means the process by 
which a novel compoimd or compoxmd composition is delivered to tiie target cell, tissue or 
organ, or placed in direct proximity of the cell, tissue, or organ. 
[0032] "Therapeutically effective amount" is the amount of novel compound or 
composition that, when administered to a subject, is effective to bring about a desired effect. 
In this case, that effect is typically an antagonistic effect that ultimately decreases the activity 
of the Na, K-ATPase. 

[0033] ^Thaimaceutically acceptable carriers" for the novel inotropic compositions may 
include sugars, starches, cellulose, powdered tragacanth; malt; gelatin; talc; excipients, oils, 
glycols, esters, agar; buffering agents, and other non-toxic compatible substances employed 
in pharmaceutical formulations. The formulations may conveniently be presented in unit 
dosage form and may be prepared by any methods well known in the pharmaceutical arts. 
Formulations suitable for oral administration may be in the form of capsules, cachets, pills, 
tablets, lozenges, powders, as granules, bolus, electuary, or a paste, as a solution or a 
suspension in an aqueous or non- aqueous Uqmd, as an oil-in-water or water-in-oil liquid 
emulsion, as an elixir or syrap, or as pastilles, each containing a predetermined amount of one 
or more of the novel inotropic compositions as an active ingredient. 
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[0034] In solid forms for oral administration, capsules, tablets, pills, powders, granules and 
the like, may by used. In the case of capsules, tablets and pills, the phamiaceutical 
compositions may also comprise bufifering agents. Compressed tablets may be prepared using 
binder, lubricant, inert diluent, preservative, or disintegrant. Liquid dosage forms for oral 
administration include pharmaceutically acceptable emulsions, microemulsions, solutions, 
suspensions, syrups and elixirs. Formulations for the topical or transdermal administration of 
a novel inotropic composition include powders, sprays, ointments, pastes, creams, lotions, 
gels, solutions, patches and inhalants, such as by aerosol. 

[0035] Formulations suitable for parenteral administration comprise a novel inotropic 
composition in combiaation with one or more pharmaceutically acceptable sterile isotonic 
aqueous or nonaqueous solutions, dispersions, suspensions or emulsions, or sterile powders 
which may be reconstituted into sterile injectable solutions or dispersions. The formulations 
may contain antioxidants, buffers, bacterostats, solutes that render the formulation isotonic 
with the blood of the intended recipient or suspending or thickening agents. Formulations 
suitable for parenteral administration may also contain adjuvants such as presCTvatives, 
wetting agents, emulsifying agents and dispersing agents, various antibacterial and antifungal 
agents, or agents that delay absorption such as aluminum monostearate and gelatin. 

Aspects 

[0036] A novel pharmacophore is disclosed. The pharmacophore allows the generation 
of basic scaffold molecules, which can then be transformed into more detailed scaffold 
molecules having inotropic properties. The pharmacophore is created to generate Na, K- 
ATPase ligands and is defined by the parameters of Table 4 and Table 5 and the spheres in 
Figure 4. The more detailed scaffolds can then be used to create potential drugs, which can 
also be tested within the pharmacophore model for blading potential. The scaffolds at least 
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satisfy the basic three point phannacophores, but may satisfy the four or five point 
phaimacophores. 

[0037] The Na, K-ATPase inhibitor scaffold compounds may be of flie formula: 




for example, wherein Rl, R2, R3 and R4 can be any organic functional group containing a 
hydrogen bond donor or a hydrogen bond acceptor and X is any element or group that allows 
the compound to retain inotropic activity. Preferably, X is a smgle element, such as N, O, S, 
or C 

[0038] The scaffold molecules are the basis functional dmgs for Na, K-ATPase 
inhibition, which can be used to prevent or treat heart conditions, particularly congestive 
heart failure and supraventricular arrhythmia. Eight compounds have been tested shown as 
inotropic agents, with three of the eight drugs have demonstrated particularly high potential. 
These drugs are shown below: 

(1) (+)-Catechin Hydrate (experimentally noted as Compound 9) 
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(2) 4'-5' Di nitroflorescein (experimentally noted as Compound 14). This compound 
demonstrates strong ATPase inhibition. 




4'.5' Dinitroflorescein 



(3) Emodin (experimentally noted as Compound 15). 




(4) Lacmoid (experimentally noted as Compound 19). This compound demonstrates strong 
ATPase inhibition. 




#19 F^NHsOrOH 



(5) R"S Dobutamine (experimentally noted as Compound 25). 
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(6) Compound 26 



(7) Compound 27 




(8) Compound 28. This compound demonstrates strong ATPase inhibition. 




The data for the tests of the eight compounds is found in Tables la and lb. 




R1 « OH or NHa 
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Table la 
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192410 


9703 






13386 


32.9 


28 


358753 










358753 


881 


9 
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Table lb 
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4.6 
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ouabain 
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12 


10 


9 


0.9 
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8692 


869 


19 


3554 


3246 




3400 


340 


28 


75572 


59710 




67641 


6764 



[0039] The drugs foimd to be efficacious in treating heart diseases will be administered to 
individuals in need thereof in pharmaceutically effective amounts for the duration of the 
disease, until the disease has abated, or as needed to prevent the disease. The novel inotropic 
compounds should also be administered in a pharmaceutically acceptable carrier. There are 
many methods and routes of administration and possible carriers, which are discussed in 
detail in the definitions section. An ordinarily skilled clinician will be able to detemiine the 
appropriate amount of the one or more novel inotropic compoimds to administer to an 
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individual. Preferably, the novel inotropic compoxmd has a wider llierapeutic index that the 
commonly used cardioglycoside dmgs, such as dioxin or ouabain. 

[0040] Finally, a general method of making and using the pharmacophore to create an 
Na, K-ATPase inhibitory compound comprises: 

(a) creating alignment between SERCAla and Na, K-ATPase, wherein SERCAla acts as the 
structural template, and the aligmnent comprises use of both dynamic programming and 
threading techniques; 

(b) transferring coordinates from a template to a model for structurally conserved regions; 

(c) generating variable regions; 

(d) refining the model through energy minimization steps; and 

(e) performing docking analysis of prospective drug candidates. 
The process may further comprise the steps of: 

(f) delineating the essential pharmacophoric elements for high binding affinity; 

(g) searching databases of known compounds using tiie restraiuts as implicated by the 
pharmacophore with allowable tolerances (e.g. ilO%); and 

(h) utilizing de novo rational drag design and computer aided molecular modeling to design 
novel compounds usiag the restraints as implicated by the pharmacophore with allowable 
tolerances (e.g. ±10%). 

[0041] As now explained in further detail, the present invention has found a way to 
utilize the recent solution of the structural coordinates of the Ei conformation of SERCAla 
(skeletal muscle sarcoplasmic reticulum/endoplasmic reticulum Ca^^-ATPase) as a structural 
template for homology three-dimensional modeling of the Na, K-ATPase. SERCA also 
occurs in the E2 conformation, and in SERCAlb, and SERCA2 forms. The Na, K-ATPase 
and other P2-type ATPases are thought to share important stmctural features with SERCAla 
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and sequence analysis of the family members established that these structures are 
evolutionarily related. 

[0042] Accordingly, by the methods of the instant invention, a model of tiie extracellular 
and transmembrane domains of the human la-subunit of the Na, K-ATPase has been 
developed, starting from the high-resolution crystal structure of SERCAla. The first and 
most important step in any modeling process is the aligmnent of the protein sequence to be 
modeled to that of the structural template. Witii a sequence identity of only -25%, aligning 
the la-subunit of the Na, K-ATPase to SERCAla by virtue of residue identity/homology is 
thwarted with complications and thus such an aligmnent cannot be confidently used as a basis 
for a tbree-dimensional stmcture. 

[0043] Therefore, in one prefCTred embodiment, the invention takes advantage of tiie 
known structural and functional similarities of the family of P2-Type ATPases and uses a 
technique that essentially ^threads" the sequence of interest to libraries of known protein 
folds from resolved structures. In this manner, the possible tertiary structure configurations 
that provide environments tiiat satisfy flie primary sequence of the la-subimits of the Na, K- 
ATPase can be investigated. Througji the independent aligmnent of the sheep la -subunit 
sequences to SERCAla, the instant invention has achieved aligmnents that have a high level 
of confidence. Areas of high local sequence similarity (indicated by shaded black and gray 
boxes. Figure 1) are well aligned as exemplified by the perfect aUgnment of the 7 residues 
(DKTTGTLT) (residues 351-357 of SEQ ID NO:l) that comprise the E1-E2 ATPases 
(aspartate) phosphorylation site that is ubiquitous to Pa-type ATPases. The TGES (residues 
181-184 of SEQ ID NO:l) sequence prior to tiie phosphorylation site and the GDGXNDXP 
(where X is any residue) (residues 701-708 of SEQ ID NO:l) motif (Figure 1, open black 
boxes) are fiirther evidence of the aligmnent. These data demonstrate that after aligning the 
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sequences without bias to residue identity, the resulting orientation effectively aligns these 
essential regions. 

[0044] Further extrapolating from the structural conservation of P2-type ATPases, the 
present invention uses computational homology modeling techniques to derive a three- 
dimensional structure for the extracellular and transmembrane domains of the la subimits of 
both human and sheep Na, K-ATPases. Thus, in another preferred embodiment of the 
invention, these three-dimensional models and specific mutations within the extracellular 
domain as they relate to ouabain affinity/binding are applied to examine the mode of 
inhibitor-receptor association. As ouabain and the other cardioglycosides bind to the 
extracellular domain of the la-subunit, five stracturally conserved regions (SCRs) gre 
modeled which collectively consist of the ten transmembrane domains and connecting 
extracellular loop regions. These SCRs are notated by sohd niunbered lines in Figure 1 and 
the numbers correlate with the helices depicted in Figure 2. Areas with extended regions of 
hydrophobicity in the primary sequence are then compared to the transmembranes predicted 
by the aligmnent, tiiose predicted by secondary structure predictors, and the transmembrane 
domains proposed as a resxilt of hydrophobicity studies. The results of the analyses are 
tabulated (Table 2) and the high level of consistency provide additional evidence to 
corroborate the alignment scheme of the invention. Furthermore, when a standard pairwise 
aligmnent is performed xising a hydrophobicity matrix, the fundamental alignment (Figure 1) 
is maintained. Consequently, supportive data firom three independent methods as well as 
structural motif conservation are obtained, all of which converge on the final sequence 
alignment of the present invention. 



17 



wo 2004/043384 



PCT/US2003/03S636 



Table 2. 
Helices 



Alignment 



Hydrophobicity Plots 



Memsat 

89-110 
123-142 
287-305 
314-338 
771-791 
800-817 
842-866 
914-930 
946-964 
978-994 



1 
2 
3 
4 
5 
6 
7 
8 
9 
10 



83-109 
126-149 
285-298 
316-337 
768-789 
799-819 
845-871 
909-934 
948-969 
981-997 



89-110 
123-140 
282-304 
311-342 
770-791 
796-818 
841-865 
909-930 
945-970 
975-996 



[0045] In another preferred embodiment of the invention, the structural coordinates of the 
helical regions of sheep la-subimit are assigned, and loops are generated that are sterically 
acceptable and that result in reasonable phi and psi angles for the spliced regions (Insightll, 
Accebys). The resulting structure is then subjected to a series of refinement steps. Initially, 
the entire protein is relaxed and theimodynamically energy inirumized to reduce initial steric 
bumps resulting from residue replacement within the secondary structural features. Next, the 
residue side chains are checked for {qjpropriate rotomer positioning using a standard rotomer 
library defined from known crystal structures. Finally, the model is energy minimized 
through a series of stq)s in which hydrogens, side chains and finally, the entire molecide are 
successively allowed freedom of movement Ribbon representations of the cioss-section and 
the extracellular side views of the final sheep model are shown in Figures 2A and 2B 
respectively. 

[0046] The resulting model can be structurally validated by the use of the wealdi of 
biological data available for the sheep la-subunit for structural validation. Site directed 
chemical labeling of the extracellular loops of the sheep Na, K-ATPase has been utilized to 
suggest structural orientation. Cysteine residues for PI 18, T309, L793, L876, and M973 
have previously been introduced to investigate membrane topology. All of the introduced 
cysteine residues could be chemically labeled indicating that these residues are exposed to the 
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extracellular surface and our model places all five residues in extracellular loop regions. 
Studies have shown that the preference of tryptophan residues and to a lesser extent lysine 
residues for membrane surfaces is a feature ubiquitous to membrane proteins. Analysis of the 
positioning of these residues in the model of the present invention shows that in fact these 
amino acids do tend to associate close to the proposed transmembrane surfaces, providing 
further vaUdation of the model. Moreover, as expected, there appears to be a propensity for 
tryptophan residues to be localized on the extracellular surface while lysine residues ^pear 
to be incorporated on the intracellular membrane. In summary, both the aligmnent and model 
of the present invention are strongly supported by experimental data. 
[0047] Jn anoflier preferred embodiment, with the information obtained &om the three 
dimensional structures of the transporter, the invention provides a pharmacophore model for 
the Na, K-ATPase, which describes the types of atoms or pharmacophoric elements and their 
geometric arrangements that are common across the set of active cardioglycosides in relation 
to the constraints imposed by the molecular binding site. The three-dimensional model of the 
pharmacophore of ravention makes it possible for the first time to examine the molecular 
mechanism of Na, K-ATPase inhibition, to screen for novel inotropic inhibitors of the Na, K- 
ATPase with higher therapeutic indices, as well as to design such inhibitors through the use 
of rational drag design techniques. With the 3D-model of the Na, K-ATPase provided by the 
present invention, the mode of cardioglycoside binding can now be elucidated. Three- 
dimensional quantitative structural-activity relationship models (QSAR) based on the 
inhibition of the sheep la-subunit as a result of incubation with cardioactive compoimds 
previously developed and the resulting comparative molecular field analysis (CoMFA) 
models provide structural implications for the binding site based on the correlation of 
changes in hgand conformation and biological activity. In combination with the model of the 
present invention, these models suggest that the length of the binding pocket is -20A, which 
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is corroborated by earlier estimates of ~19A. Mutational studies suggest that the ligand 
binding site is situated between H1-H2 (Ribbon No. 1, Figure 2A) and H5-H6 (Ribbon No, 3, 
Figure 2A), and H1-H2 and H9-H10 (Ribbon No. 5, Figure 2A) corresponding to distances of 
20A and 22A respectively. CoUectively these data suggest that the binding site consists of a 
--20A groove comprised of the extracellular loops connecting H1-H2, H5-H6, and H9-H10. 
[0048] Thus, in another preferred embodiment, the present invention provides a three- 
dimensional model of the Na, K-ATPase cardiotonic binding site. In a fiirther preferred 
embodiment of the invention, a pharmacophore model is provided, which precisely defines 
the mechanism of interaction between the Na, K-ATPase and an inhibitory molecule. This 
model allows for the de novo design of drug candidates based on the Na, K-ATPase receptor 
structure. The pharmacophore allows for the development of inotropic agents that are 
structurally unique from those chemical classes known and exannned to date, thus making it 
possible to develop inotropic agents with higher (also called wider) therapeutic indices. 
[0049] The systematic docking of 16 cardioglycosides inside the putative Na, K-ATPase 
binding pocket provided by the instant invention was achieved by using ) GOLD (Genetic 
Optimization for Ligand Dockmg) Jones et al JMolBiol 254, 43-53 (1995). However, any 
accurate docking program is contemplated, including conformational sampling-based 
methods (e.g., free energy perturbation, linear interaction energy approximations), 
knowledge-based ''potential of mean force" methods (Pmf, Dragscore), force-field methods 
(GOLD, Dock, AutoDock), Apropros, CombiBUILD, Escher, FlexiDock, FLOG, GRAMM, 
HEX, HotDock, LIGIN, PUZZLE, STALK, and empirical free-energy scoring fimctions 
(Ludi, Chemscore, Score, Fresno, FlexX, Pip). In addition, the docking and scoring methods 
of US Provisional AppUcation Number 60/477,714, entitled •'Novel Knowledge-based 
Docking and Scoring Algorithm for Virtual Drug Screening (Welsh, et al.), filed June 10, 
2003 is incorporated by reference herein. 
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[0050] GOLD is a ligand-docking program, which uses a genetic algorithm (GA) to 
explore ligand confonnation and satisfy ligand-binding requir^ents 
(http://www.ccdc.cam.ac.uk/prods/gold). One advantage over other docking systems is that 
GOLD allows not only for ligand flexibility, but also for limited flexibility of the binding 
pocket creating a more realistic environment for inhibitor association. Further, the stochastic 
nature of GOLD ensures that the search space is well explored and local extremes are less 
influential as compared to deterministic algorithms. Thus, it was possible to select a subset 
of inhibitors consisting of all molecules 15-30A in length and including ouabain, digoxin and 
digitoxin. These molecules best fit the binding site of the invention, and as a result are 
imlikely to have more than one mode of association. The combination of experimental and 
computational results suggests that the ligand binding domain is '-20A which effectively 
defines a docking sphere of a --lOA radius. A centroid of the sphere was determined as the 
intersection of the vectors defining H1-H2 with H5-H6, H1-H2 with H9-H10 and H5-H6 with 
H9-H10. Each inhibitor was docked 25 independent times. 

[0051] To explore more thoroughly the protein-inhibitor interaction, the radii of the 
docking spheres were varied (with increasing radii lOA, 12A, and 15A). Similar results were 
obtained with the two smaller docking spheres while the excessively large volume of the 15 A 
docking sphere resulted in unrealistic interactions in which the molecules docked 
perpendicular to the molecular surface of the Na, K-ATPase. Using the results generated 
with the 12A-docldng sphere, the consensus mode of binding for each inhibitor from the 
family of 25 dockings was detemiined. As all the binding modes were rank ordered using the 
GOLD fitness function, which measures steric and electrostatic complementarity between 
each ligand conformation and the receptor, the molecule corresponding to the highest scoring 
confomiation of each binding consensus was selected. As an example, the aligmnent of the 
25 conformations of ouabain (Figure 3 A) and the overall alignment of all 16 inhibitors 
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representing the consensus conformation (Figure 3B) is provided. It is notable that while 
GOLD is a stochastic algorithm, the independent dockings converged to a single consensus 
both within dockings of an independent compound as well as across dockings between 
compounds. In a manner consistent with experimental findings, the cardiotonic steroid 
moiety and the unsaturated lactone ring are oriented toward H1-H2, while the carbohydrates 
are directed toward H9-H10. Compounds with more than one sugar see these moieties wrap 
around H7-H8. All compounds completely occupy the potential binding pocket which is 
situated over the channel (depicted by a fast Connolly channel surface) thus effectively 
inhibiting the function of the ATPase by blocking the ion channel (Figure 3C). 
[0052] The model of the binding site allows for two primary observations: First, that the 
consensus binding orientation is highly conserved. With no outliers, the dockmg modeling 
results in a single binding mode and suggests that only the cardiotonic steroid moiety, lactone 
ring and alpha sugar contribute to binding. Second, while the overall domain of binding is 
very consistent, the individual orientation of the steroid backbones of the inhibitors is not 
uniform indicating that this moiety is acting as a scaffold which provides appropriate spacing 
for the functional groups, rather than as a phannacophoric element essential to binding. 
Furthermore, the number of sugars has no effect on the binding consensus. For example, 
digitoxin with two or three sugars maintains an identical binding orientation (Figure 3D). 
This observation agrees not only with previous in vitro studies describing the orientation of 
the sugars of bound inhibitors, but furthermore corresponds with our own observations based 
on our CoMFA models. 

[0053] One factor that influences the Gibb's free energy (AG) of binding of a small 
molecule to a receptor is the formation of hydrogen bonds between the compound and the 
receptor. One of the most widely used definitions to determine hydrogen bonding potential 
for secondary structure determination is the Database of Secondary Stmcture in Proteins or 
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DSSP. Systematic measurement of the distance between all potential hydrogen bond donors 
and/or acceptors of flie ligands and the residues comprising the binding pocket of the 
invention, the potential residue interactions of each compound within the binding site were 
m^ped. Of the 16 molecules docked to the la-subunit, 15 had the potential to fomi 
hydrogen bonds with the side chain or backbone nitrogen of residue D121 . Biological data 
has suggested that D121 residue is essential for inhibitor binding as a single point mutation 
(D121N) at this location decreases sensitivity to ouabain by 1000 fold. Four additional 
residues also appear to be involved in molecular inhibition: the side chains of Ql 1 1 and 
N122, and the backbone nitrogens of E908 and M973. The data provided by the methods and 
models of the present invention show that 15 of the 16 molecules docked could potentially 
form hydrogen bonds with Ql 1 1, and mutants of tihie la-subunit of the Na, K-ATPase with R 
replacing residue Ql 1 1, or R and D replacing residues Ql 1 1 and N122 respectively, are 8.3 
fold and 1 150 fold less sensitive to ouabain than the wild-type sheep catal3^c subunit. In 
addition, all 16 molecules were able to forai hydrogen bonds with the backbone nitrogen of 
E908, and 12 molecules were shown to have the potential to interact with M973. 
[0054] The effects of these mutations on the binding of the cardioglycosides can further 
be predicted by modeling the la-subunit of the Na, K-ATPase incorporating either the 
Ql 1 IR or D121N mutations. In addition, the mutations were introduced in the models of the 
human and rat la-subunits of the Na, K-ATPase whose sequence differs from that of the 
sheep by 23 residues (or 2%) and 35 residues (3.5%) respectively. The human la-subunit 
remains sensitive to ouabain while the rat subunit is insensitive to ouabain inhibition. The 
docking study with ouabain was repeated for both the human and rat la-subunits, and the 
mutant structures maintaioing identical parameters as previously. Ouabain binding to the 
human la-subunit of the Na, K-ATPase is identical to that observed in the original docking 
study. The binding affinity of ouabain for both the human and the sheep la-subunit of the 
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Na, K-ATPase is analogous, thus inclusion of the human model acts as a control for the 
following mutational studies. The incorporation of either single point mutation (Ql 1 IR or 
D121N) results in the reorientation of the lactone moiety from H1-H2 toward H3-H4. 
Furthermore, the rat (which contains both Ql 1 IR and N122D point mutations) forces 
ouabain to abandon the consensus orientation observed with the human and sheep subunits. 
[0055] The mutational effects with respect to ouabain sensitivity, described above, can be 
further explored computationally as a decrease in receptor sensitivity to an inhibitor often 
reflects diminished binding affinity. As such, binding energies (AEbmding) associated with 
various ouabain-enzyme interactions resulting from the docking studies were calculated and 
then compared with available inhibition data. In order to quantify the reorientation of ouabain 
with respect to the affinity decreases observed experimentally, the relative change in binding 
energy (AEbinding) were calculated. Values of AEbmding were computed as the difference in 
calculated energy between the ouabain-receptor complex (Ecompiex) and the simi of the 
energies calculated separately for the free ligand (Eiigand) and the free receptor (Exceptor). 
[0056] The molecular mechanics calculations included bonded and nonbonded (vdW and 
Coulombic) energy terms and the electrostatic component of solvation. The AEbmding data 
(Table 3) shows a strong statistical correlation with the experimentally derived inhibitory 
effects of the mutants (Pearson correlation r = 0.8 1 ; Spearman correlation r = 0.88). The 
calculated values of AEbinding for the wild-type sheep and human al-subunits are consistent 
with their nearly equivsdent ouabain affinity. Furthermore, introduction of the Ql 1 IR 
mutation which is known to lower ouabain's affinity for the sheep enzyme by -10-fold 
resulted in a corresponding decrease in the calculated binding affinity. Furthermore, the rat 
al-subunit, or introduction of the D121N into the sheep oil model, yielded AEbinding values 
consistent with the dramatic loss (-1000 fold) loss in ouabain inhibition associated with 
either of these enzymes. 
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[0057] As the results of our molecular modeling studies have impUcated two additional 
residues as important components of the digitalis binding site, it was sought to further 
investigate the role of these amino-acid residues with respect to ouabain binding. While in 
vitro mutation studies are beyond the scope of tiie present study, these mutations were 
simulated in silico to determine the influence of these residues on the affinity of ouabain. 
Mutant sheep al-subunits were constructed incorporating either an E908L or M973A 
mutation and docking analyses with ouabain were performed. The consensus orientation of 
ouabain selected for analysis was obtained as described above. 

[0058] Interestingly, the consensus orientation observed from the docking ensemble 
differed from the consensus orientation of ouabain bound to the wild-type sheep enzyme 
suggesting that these residues each indq)endently play a role in ouabain binding. Values of 
AEbinding were calculated (Table 4) and compared with the AEbmding values calculated for the 
wild-type and mutant sheep-, rat- and human enzymes. The results indicate that both 
mutants could have a large impact on ouabain binding affinity. 
Table 3. Calculated changes in Binding Energies (AEbindmg) 



Species/Mutant* 


Relative Inhibition 


AAEbindine (Iccal/mol) 


Sheep 


1 


-204 


Rat 


1000 


-158 


Human 


1 


-199 


D121N 


1000 


-107 


QlllR 


10 


-167 


E908L 


ukn 


-160 


M973A 


ukn 


-131 



^Mutation of sheep structure. 



[0059] In another preferred embodiment, the invention provides a method of identifying 
novel ligands of the Na, K-ATPase. Identification of novel ligands can be achieved by 
methods known in the art, such as the screening of libraries and databases for chemicals 
including proteins, peptides, nucleotides, and small molecules that correspond to the 
stractural criteria defined by the pharmacophore of the invention. The pharmacophore of the 
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present invention contains the essential hydrophobic and hydrophilic regions for interaction 
(Figure 4). 

[0060] There are essentially three regions of the binding pocket that could be essential for 
binding: a hydrogen bond donor/acceptor corresponding to the potential hydrogen bond 
donor/acceptor of Ql 1 1, an intervening hydrophobic area, and a hydrogen bond donor 
corresponding to the potential hydrogen bond acceptor D121 (represented respectively by 
spheres 1,2 and 3; Figure 4). The two hydrogen bonding groups are 9.95 A and 7.63 A from 
the hydrophobic core, resulting in a total distance of 17.10 A separating the two hydrogen 
bondmg (Figure 4A). In addition to the three pharmacophoric elements discussed above, the 
3D coordinates of E908 (sphere 4; Figure 4) and M973 (sphere 5; Figure 4) provide 
positioning for auxiliary potential hydrogen bond donor/acceptor groups. The rigor of the 
database search criteria can thus be increased in a stepwise fashion by introducing four and 
five point pharmacophores (Table 4). 



Table 4. 



Pharmacophore 


Sphere 1 


Sphere 2 


Sphere 3 


Sphere 4 


Spheres 


Three Point 


HBD 


HYD 


HBD 






HBA 


HYD 


HBD 






Four Point 


HBD 


HYD 


HBD 


HBD 






HBA 


HYD 


HBD 


HBA 






HBD 


HYD 


HBD 




HBD 




HBA 


HYD 


HBD 




HBA 




HBD 


HYD 


HBD 


HBA 






HBA 


HYD 


HBD 


HBD 






HBD 


HYD 


HBD 




HBA 




HBA 


HYD 


HBD 




HBD 


Five Point 


HBD 


HYD 


HBD 


HBD 


HBD 




HBD 


HYD 


HBD 


HBD 


HBA 




HBD 


HYD 


HBD 


HBA 


HBD 




HBD 


HYD 


HBD 


HBA 


HBA 




HBA 


HYD 


HBD 


HBD 


HBD 




HBA 


HYD 


HBD 


HBD 


HBA 




HBA 


HYD 


HBD 


HBA 


HBD 




HBA 


HYD 


HBD 


HBA 


HBA 



*HBA represents Hydrogen Bond Acceptor; HBD represents Hydrogen Bond Donor; HYD represents 
Hydrophobic group 
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Distances between points are located in Table 5. 



Table 5. 



Pharmacophoric 
Elements 



2 



3 



4 



5 



1 
2 
3 
4 



15 ±0.75 



8 ±0.40 18 ±0.90 22 ±1.10 

9 ±0.45 6 ±0.30 9 ±0.45 

- 12 ±0.60 14 ±0.70 
7 ± 0.35 



*A11 measurements are in A. 



[0061] However, compounds do not have to satisfy all criteria to be active in inhibiting 
the la subunit of the ATPase. The combination of the initial three pharmacophoric elements 
(spheres 1, 2, and 3) collectively represent the minimum 3D criteria necessary to search 
existing databases for novel potential scaffolds. To validate the basic three-point 
pharmacophores of the invention, NCI and Maybridge three-dimensional databases were 
searched for molecules that satisfy the pharmacophoric parameters. The pharmacophores of 
the invention allowed the identification of existing cardioglycosides that were independent of 
those used for pharmacophore development in the databases utilizing the pharmacophore. 
Therefore, the three-dimensional selection criteria developed by the present invention mirror 
the physical elements required for inhibitor binding and can be therefore used for the design 
of new therapeutic agents. In addition to the identification of cardioglycosides, the database 
search identified several new potential scaffolds of interest. Based on the models and 
methods of the invention, it was shown that the cardiosteroid backbone of the inhibitors are 
acting merely as hydrophobic spacers, and therefore substitution of scaffolds that maintain 
the three dimensional orientation of the pharmacophoric elements should succeed in 
inhibiting the Na, K-ATPase. 

[0062] Thus, in yet another preferred embodiment of the invention, a method of 
designing ligands that will inhibit the Na, K-ATPase is provided. The structural data 
provided by the models and methods of the present invention make it possible to elucidate the 
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mechanism of cardioglycoside inhibition of the la-subunit of the Na, K-ATPase. This in 
turn allows the design ofde novo inhibitors of the la-subunit of the Na, K-ATPase. Several 
computer-based schemes using molecular modeling software may be employed to build 
potent and selective inhibitors from these scaffolds. Candidate scaffolds can be placed in the 
binding site of the target molecule, and a Ugand can be built by successively bonding other 
fragments to it using methods known in the art. Thus, the Na, K-ATPase model and 
pharmacophores of the present invention will be used to guide the * growth' of a potential 
inhibitor. Public and conamercial chemical Ubraries (ACD, Maybridge, NCI, and WDI) may 
be searched for new molecides or scaffolds that satisfy the query's pharmacophore withia 
specified tolerances. Further, computational tools will be employed to match small 
molecules with each other or with receptor topology rapidly based on similarity or 
complementarity iu shape and electrostatic properties. 

[0063] Promising candidate inhibitors emerguig from the computer-aided design pipeline 
can be evaluated in terms of synthetic feasibility and overall acceptabiUty. Following 
chemical synthesis, binding and inhibition assays known in the art can be used to evaluate the 
in vitro effects of the potential inhibitors. 

Experimental Details 

[0064] The following detailed examples are provided to illustrate the claimed invention 
better and are not to be interpreted as limiting the scope of the invention. To the extent that 
specific materials are mentioned, it is merely for purposes of illustration and is not intended 
to limit the invention. One skilled in the art may develop equivalent means or reactants 
without the exercise of inventive capacity and without departing from the scope of the 
invention. 
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[0065] The sequences of the al-subunits of the human, sheep, and rat Na^ K^-ATPases 
were obtained from GenBank (Accessions nos. P05023, P04074, P06685 respectively). The 
SERCAla crystal structure (lEUL) was obtained from RCSB-PDB and the sequence was 
extracted using Insight n (Accehys, San Diego, CA), The sequence alignment between 
SERCAla and the al-subunit of the sheep Na^ K^-ATPase was deduced from matrix based 
(PAM250 and Blosum35) approaches (www.ebi.ac.uk/clustalwA and by threading the sheep 
sequence to the structural features of SERCAla using THREADER 

(http:/^ioinf cs.ucl.ac.uk) from University College London and the resulting aUgmnent was 
ftirther manually refined. The amino acid identity resulting from this alignment is 25,6% 
and all local areas of structural and functional importance are well conserved. The 
extracellular gaps resulting &om the alignment are localized to loop regions (H1H2, and 
H7H8) and are thus easily accommodated. 

Homology Modeling of the Na. K-ATPase 

[0066] All molecular modeling operations were carried out using Insightn and related 
modules (Accelrys, Ihc,, San Diego, CA). The three-dimensional coordinates of the 
extracellular and transmembrane domains of the al subunit of sheep Na^ JC^-ATPase were 
assigned. Undefined loops and gaps in the aligned sequences were accormnodated by 
random generation of intervening sequences while conserving the integrity of backbone 
orientations in stracturally defined, flanking residues. Inclusion of the randomly modeled 
regions required splice site relaxation in order to allow the protein backbone to adopt 
appropriate psi/phi angles. This operation was completed by relaxing the backbone atoms of 
those residues involved in spUce sites and subjecting the molecule to a localized energy 
minimization procedure. 
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[0067] After generation of the model, the side chain residues were checked for acceptable 
rotomer positioning using the internal rotomer library and steric clashes were resolved by 
brief energy minimization of the entire molecule (250 iterations, steepest descent). Finally, 
the model was subjected to a fall energy minimization through a series of indepradent steps 
in which hydrogen atoms, side chains, backbone atoms and, ultimately, the entire molecule 
were successively allowed freedom of movement. Similarly, the human and rat sequences 
were aligned to the sheep sequence, and then modeled and refined as described above. Point 
mutations as described herein were incorporated into the faUy refined sheep structm*e, and the 
resulting mutants were again energy minimized through the steps discussed above. This 
procedure implemented the CVFF force field with a distance dependent dielectric fanction 
(8=Sor, with So = 4) until the convergence criterion of 0,04184 kJ/mol (0.01 kcal/mol) change 
in energy between successive iterations was achieved. The dielectric constant = 4 was 
chosen as it corresponds to that of water in a fast-moving dynamic field. 

Development of the Pharmacophore Model of the Na, K-ATPase 

[0068] The 3D structure of each cardiotonic compound was built using the molecular 
firagment library of Sybyl6.6 (Tripos, Inc., St. Louis, MO) and based on the structure of 
digoxin as determined by X-ray crystallography.(ref) The conformation of the P- and y-sugars 
was obtained by global energy miiiiniization detemiinations. As previously described, the 
final structural models were energy minin^zed using the Tripos force field with the 
Gasteiger-Marsili method for assigning partial atomic charges.(reQ The GOLD operator 
usage probabilities for crossover, mutation, and migration were set to 95%, 95%, and 10%, 
respectively. The maximum distance between hydrogen bond donors and fitting points was 
set to 5A, and the nonbonded cutoff for the van der Waals energies was adjusted to the 
manufacturer recommended settings (lOk/y = well depth of the van der Waals energy of the 
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atom pair ij). The ligand orientation accepted was the top scoring consensus orientation of 
25 independent genetic algorithm (GA) runs, each with a maximum number of 1000 GA 
operations perfomied on a single population of 100 individiials. The selection pressure was 
1.1. The centroid of the binding sphere was determined as the intersection of the vectors 
connecting H1H2 witii H5H6, H1H2 with H9H10 and H5H6 with H9H10. 

[0069] Binding energy calculations were performed as follows: AEbmdmg defined as the 
difference between the potential energy of the ouabain-receptor complex (Ecompiex) and the 
sum of potential energies of the ouabain (Eiigand) and receptor (Ereceptor). 

AEbinding ~ Ecomplex ~ (Eligand Ereceptor) 

Energy values represent the niinimum potential energy resulting firom the optimization 
procedure described above. Ouabain was allowed full freedom of movement througihout the 
process. Various orientations of ouabain with respect to the ATPase were obtained by 
combining the consensus orientations of ouabain as determined by the GrOLD docking, with 
the appropriate homology modeled enzyme. Specifically, the bonded (Stretch, torsion etc.) 
and nonbonded (van der Waals and Columbic terms) energy terms were calculated by 
Discover (Jhsightll) whilst the linear Poisson-Boltzmaim equation was solved using DelPhi to 
obtain the electrostatic contributions of the solvation energy (reaction field energy). 

Design of Novel Na. K-ATPa< se TnTiihitnrs 

[0070] Candidate inhibitors identified from the receptor-based pharmacophores provided 
by the present invention will be further screened in terms of bioavailability and toxicological 
profile. Bioavailability is a measure of the value of the hydrophobic/hydrophilic balance, as 
measured by log?, will be predicted using the Clog? program (Daylight Chemical 
Information Systems, Santa Fe, NM). Values of the aqueous solubility and cell membrane 
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permeability will also be predicted using available molecular models in the laboratory as 
additional screens to select only the most promising candidate inhibitors. The toxicological 
profile of potential inhibitors will be determined by calculating relevant endpoints including 
carcinogenicity, Ames mutagenicity, developmental toxicity potential, skin irritancy, ocular 
irritancy, and aerobic biodegradability will be predicted for the candidate inhibitors using the 
QSAR-based TOPKAT program (Accehys, hic, San Diego, CA). 

[0071] Several computer-based schemes may be employed to build potent and selective 
inhibitors from these scaffolds* First, a method known as "molecule growth" may be 
employed to place a *seed' fragment or scaffold in the binding site of the target molecule, and 
to build a ligand by successively bonding other fragments to it using methods such as Small 
Molecule Growth (SmoG) (DeWhitte and Shakhnovich, J. Am. Chem. Soc. 118: 11733- 
1 1744, 1996)) or GrowMol (Bohacek and McMartin, J. Am. Chem. Soc. 1 16: 5560-5571, 
1994). The Na, K-ATPase homology model of the present invention will be used to guide the 
'growth' of a potential inhibitor. 

[0072] Second, pubUc and commercial chemical Ubraries (ACD, Maybridge, NCI, and 
WDI) may be searched using a fragment-based development method commonly known as 
"sub-structure searching." This method uses a 3D pharmacophore as a query to search a 
chemical library for new molecules (or scaffolds) that satisfy the query's pharmacophore 
within specified tolerances. Third, a computational tool will be employed as a ray-tracing 
scheme to rapidly match small molecules with each other (or even with receptor pockets) 
based on similarity (or complementarity) in shape and electrostatic properties. It is a 
powerful and efficient method for dmg discovery and particularly for the discovery of new 
scaffolds ttiat share the same properties but not necessarily the same structural formula as the 
query molecule. 
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[0073] The design strategies used for chemical synthesis will depend on the structures of 
the proposed inhibitors. Promising candidate mhibitors CTierging from the computer-aided 
design pipeline will be evaluated in terms of synthetic feasibility and overall acceptability. 
Following chemical synthesis, binding and inhibition assays will be used to evaluate the in 
vitro effects of the potential inhibitors. Assays will be conducted as describe previously 
(Biochemistry 2002, 41, 1137-1 148). Briefly Na, K-ATPase will be purified from the outer 
medulla of frozen lamb kidney. Competitive binding assays will be performed in the 
presence and absence of varying concentrations of the potential inhibitors. Inhibition studies 
will be carried out with and witiiout varying concentrations of the potential inhibitors, and the 
NADPH oxidation rate in the presence and absence of the inhibitors will be determined. 
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